Critical property of spin-glass transition in a 
bond-disordered classical antiferromagnetic 
(N : Heisenberg model with a biquadratic interaction 



C/3 



o 



> 



o 



X 



Hiroshi Shinaoka^, Yusuke Tomita^ and Yukitoshi Motome^ 

^Nanosystem Research Institute, AIST, Tsukuba 305-8568, Japan 

^Institute for Solid State Physics, University of Tokyo, Kashiwanoha, Kashiwa, Chiba, 

277-8581, Japan 

■^Department of Apphed Physics, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 

113-8656, Japan 

E-maih h.shinaoka@aist.go.jp 



2 ' Abstract. Motivated by puzzhng spin-glass behaviors observed in many pyrochlore-based 

magnets, effects of magnetoelastic coupling to local lattice distortions were recently studied by 
the authors for a bond-disordered antiferromagnet on a pyrochlore lattice [Phys. Rev. Lett. 
107, 047204 (2011)]. Here, we extend the analyses with focusing on the critical property of the 
spin-glass transition which occurs concomitantly with a nematic transition. Finite-size scaling 
analyses are performed up to a larger system size with 8192 spins to estimate the transition 
temperature and critical exponents. The exponents are compared with those in the absence 
of the magnetoelastic coupling and with those for the canonical spin-glass systems. We also 
discuss the temperature dependence of the specific heat in comparison with that in canonical 
spin-glass systems as well as an experimental result. 



1. Introduction 

Spin glass (SG), in which spins are frozen in random directions, is a low-temperature(T) phase 
widely observed in geometrically frustrated magnets. Recently, various puzzling behaviors have 
been observed experimentally for the SG transition in several pyrochlore-based magnets. One is 
the robustness of the SG transition temperature Tf. For example, pyrochlore oxides -R2M02O7 
show almost constant Tf ~ 20K for R=Y, Dy, and Tb despite largely different Curie- Weiss 
C^ \ temperatures [H El El IH E]- Furthermore, in Y2M02O7, Tf is practically unchanged against 

random substitution of Y^+ by La^^" [S]. Similar behaviors are observed in another pyrochlore- 
based magnets, e.g., spinel oxides i?Cr204. In ZnCr204, random substitution of Zn^"*" by Cd^'^ 
destroys a spin-lattice ordering and induces a SG state with fairly flat Tf [7]. Another puzzling 
aspect is that Tf is much higher than a theoretically expected value for a moderate strength of 
disorder [HI [9l [TO]. These behaviors suggest that Tf is determined neither by a typical energy 
scale of exchange interactions nor the strength of disorder; some important factor is missing in 
the previous SG theories to understand the SG behaviors in these frustrated magnets. 

Stimulated by these puzzles, the authors recently developed a SG theory by including a 
magnetoelastic coupling to local lattice distortions in a bond-disordered pyrochlore Heisenberg 
antiferromagnet [11]. Using classical Monte Carlo (MC) simulations, we revealed that Tf is 



largely enhanced by the spin-lattice coupling. Moreover, we showed that Tf is set by the energy 
scale of the spin-lattice coupling, and becomes almost independent of the disorder strength of 
the exchange interaction in a wide range of parameters. These results explain well the puzzling 
aspects of the SG transition observed in the pyrochlore-based magnets. 

An interesting question is whether the critical property of the SG transition is modified by the 
magnetoelastic coupling. In the previous study, using finite-size scaling analyses, we showed that 
the MC data are compatible with a second-order transition, and obtained the critical exponents 
7 = 1.2(5) and v = 0.67(16) [TT]. These critical exponents agree with the previously-reported 
ones in the absence of the spin-lattice coupling within the error bars [51 El [TO] • The values are 
also consistent with recent numerical estimates for the canonical SG \\.'1\ I13j . These results 
suggest that the magnetoelastic coupling as well as the geometrical frustration does not modify 
the critical property. The analyses, however, were limited to small system sizes, i.e., up to 2000 
spins in Ref . [11] . It is well known that a careful consideration on finite-size effects is necessary in 
estimating the critical exponents for the canonical SG. Here, we report more elaborate analysis 
up to a larger size with 8192 spins, and discuss the critical property in more detail. 

2. Model and method 

We consider the following antiferromagnetic Heisenberg model with a biquadratic interaction 
introduced in Ref. [TT] : 

% = 2_^ JijSi ■ Sj — bij[Si ■ Sj) , (1) 

{id) 

where Si denotes a classical Heisenberg spin at site i and the sum runs over the nearest neighbor 
bonds of the pyrochlore lattice. The model is derived from a spin-Peierls type model by 
integrating out local lattice distortions. Here, Jij is a bond-disordered exchange interaction, 
randomly chosen from [J — A, J + A] with < A < J. Effects of the magnetoelastic coupling 
are taken into account by the biquadratic interaction bij] we measure its strength by 5 (= bij/jf-). 

The previous MC study showed that the model exhibits a peculiar multicritical phase diagram 
between a nematic phase and SG phase [TT]. In the weakly disordered region A < 6, there are 
successive two transitions, a nematic transition at Tc ~ 6 and a SG transition at Tf < Tc. Tf 
grows as A increases, and appears to merge into Tc at A ~ b. For larger A, Tf (= Tc) becomes 
nearly independent of A. In the following, we focus on the critical property of the concomitant 
transition in the region of A > 6. Specifically, we take 5 = 0.2 as in the previous work, and 
perform MC simulations at A = 0.8. 

In the following calculations, we consider periodic systems of cubic geometry with L^ unit 
cells with totally Ng = 16L^ spins. The largest system employed in the following simulations 
is L = 8 with 8192 spins. In the MC sampling, we adopt the exchange MC method [U] and 
the overrelaxation update [TO] in addition to the standard single-spin-flip update. To clarify 
the critical behaviors of the SG and nematic transitions, we calculate the SG susceptibility 
XsG) nematic susceptibility XQi and specific heat C. The spin-glass susceptibility XSG is 
given by Ngq'^j^, where g^^ is the Edwards- Anderson order parameter defined by q|^ = 

Ni^{{J2^,u=x,y,ziT,^=iSfl,S^:,f)T)A [IS]. Here (•••)t denotes a thermal average and (• • •)a a 
random average over the interaction sets; a and /3 denote two independent replicas of the system 
with the same interaction set. The nematic susceptibility XQ is given by XQ = ^sQ^/T, where 
Q2 is the spin collinearity defined by Q^ = 2N~'^{{Y,ij{{Si-Sjf -1/3})t)a- The specific heat C 
is calculated by C = N^^T^^ {{T-i'^)'r — {'H)'^)a- In the following MC simulations, we interpolate 
the calculated thermodynamic quantities using the histogram reweighting technique [T7]. All 
the data shown below are averaged over a number of interaction sets varying from 256 to 2000 
depending on L. 
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Figure 1. Calculated thermodynamic quantities at 6 = 0.2 and A = 0.8: (a) the spin-glass 
susceptibility XSG) (b) the nematic susceptibility xq, and (c) the specific heat C. The system-size 
dependence of the peak value of C is shown by a log-log plot in the inset of (c) . 



3. Results 

Figures 1 shows the MC results at 6 = 0.2 and A = 0.8. To examine the critical property of 
the concomitant transition of the nematic and SG, we perform finite-size scaling analyses for 
XSG and xq separately. Here we assume x = L'^''^ f{L}'^t), where t = {T — Tt)/Tt and (x, Tt) 
is either (xsG) ^f) or (xq, Tc). v and 7 are the critical exponents for the correlation length 
and X) respectively. As demonstrated in Figs. 2(a) and 2(b), we successfully obtained scaling 
collapses of xsG and XQ for 4 < L < 8: From the scaling of xsG [Figs. 2(a)], we obtained 
Tf = 0.248(2), 7 = 1.5(1), and v = 0.80(2), and from that of XQ [Figs. 2(b)], T^ = 0.249(1), 
7 = 1.6(1), and u = 0.78(2). The results are summarized in Table 1 with the estimates by using 
a set of data for smaller system sizes 2 < L < 4. Note that we observed no significant system- 
size dependence in the scaling results: The estimates for different range of L coincide with each 
other within the error bars. The results that the two transitions have the identical exponents 
and transition temperature confirm the conclusions in the previous study [11]: Two transitions 
occur concomitantly in a second-order fashion, and the critical exponents are indistinguishable 
from those for 6 = [3 El HO] and the canonical SG [T2|[T3]. 

Figure 1(c) shows the T dependence of C for various system sizes. C exhibits a single broad 
peak at r ~ Tf(= Tc), in contrast to the first-order nature of the nematic transition at A = 0. 
The peak value Cpcak shows a very weak L dependence, and its growth gets slower as L increases, 
as shown in the inset of Fig. 1(c). When we fit the data by Cpeak ex L", the exponent a decreases 
as L increases [a = 0.076(1) for 2 < L < 4, and a = 0.046(2) for 5 < L < 8]. Alternatively, the 
data can be well fitted by assuming C~ j, = aL^ + C~^{qo) with p ~ — 1 (not shown). These 
results suggest that C is non-singular in the thermodynamic limit. The broad peak behavior is 
superficially similar to that observed in the canonical SG [18], but the peak is located near Tf 
in the present case: the peak position is estimated as Tpeak = 0.248(5) for L = 8. This is in 
contrast to the case of the canonical SG in which Tpcak exceeds Tf typically by 20% [18]. The 
broad peak at Tf will be of characteristic of the SG transition concomitant with the nematic 
transition in the present system. Note that the result is consistent with the magnetic specific 
heat measured for Y2M02O7, in which a broad peak is observed near Tf ~ 20K [19] . 



4. Summary 

In this paper, we have examined the critical property of a concomitant spin-glass and nematic 
transition in a bond-disordered pyrochlore antiferromagnet with a biquadratic interaction, which 
was recently found by the authors [11]. We have extended Monte Carlo simulations and finite- 
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Figure 2. Scaling 
collapses of (a) XSG 
and (b) XQ for the 
data in the range of 
0.225 < T < 0.275. 



System size 


XSG 


XQ 




Tf 


7 


V 


Tc 


7 


u 


L = 2,3,4 

L = 4,5,6,8 


0.248(2) 
0.248(2) 


1.57(7) 
1.5(1) 


0.79(2) 
0.80(2) 


0.242(2) 
0.249(1) 


2.0(3) 
1.6(1) 


0.87(7) 
0.78(2) 



Table 1. Estimates of the transition temperatures and critical exponents. The system sizes 
employed in the finite-size scaling analyses are L = 2, 3, 4 and L = 4, 5, 6, 8, respectively. 



size scaling analyses up to a larger system size with 8192 spins. As a result, we have confirmed 
the previous report that the concomitant transition is second order, and its critical property 
does not contradict those in the absence of the spin-lattice coupling and for the canonical spin 
glass. Moreover, we have discussed the temperature dependence of the specific heat. 
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